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ABSTRACT 

We report on the measurement at 820- and 1400-MHz of orbital modulation of the diffractive scintillation 
timescale from pulsar A in the double-pulsar system J0737-3039 using the Green Bank Telescope. Fits to this 
modulation determine the systemic velocity in the plane of the sky to be Viss ~140.9±6.2kms _1 . The parallel 
and perpendicular components of this velocity with respect to the line of nodes of the pulsar's orbit are Vpi ane 
~96.0±3.7kms _1 and V P ei P ~103.1±7.7kms _I respectively. The large Vp erp implies that pulsar B was born 
with a kick speed of >100kms _1 . Future VLBA determination of the angular proper motion in conjunction 
with improved Viss measurements should provide a precise distance to the system. Using high-precision timing 
data and the Viss model, we estimate a best-fit orbital inclination of i = 88. °7 ± 0.°9. 
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1, INTRODUCTION 

For over twenty years observers have known that mea- 
surements of the decorrelation bandwidth Avj and scintil- 
lation timescale Atd of pulsars undergoing strong Diffrac- 
tive Interstellar Scintillation (DISS) ca n be used to e stimate 
their velocity in the plan e of the sky (Lvne & Smith 1982). 
ICordes & Rickettl (^998) examined DISS-derived pulsar ve- 
locities in detail and found that the measurements depend 
heavily on the observing frequency, the direction and distance 
to the pulsar, and the detailed distribution of the interstel- 
lar material causing the scintillation. This last point makes 
measurements of pulsar velocities particularly difficult since 
different models for the electron distribution and its irregu- 
larities can cause differences in the "measured" scintillation 
velocities (Viss) by factors of a few. However, for binary 
pulsars in compact orbits (orbital periods P 01 -b < 1 day), the 
pulsar-timing-derived orbital velocities can be used to cali- 
brate Viss measurements and remove many model-dependent 
and/or systematic effects. Unfortunately, suitable binary pul- 
sars are rare and successful measurements of this kind have 
only b een made for two pulsars: PSR B06 55+64 by | Lyne 
(1984) and PSR Jl 141-6545 bv lOrd. Bailes. & van S traten 
(2002, hereafter OBvS). 

PSR s J0737-3039A & B (hereafter A and B;|Bureav et al. 
120031 Evne et al J 120 04) comprise a fantastic double-pulsar 
binary (hereafter 0737) consisting of the 22.7-ms pulsar A 
and the 2.77-s pulsar B. It is nearby (^0.6 kpc), mildly ec- 
centric (e ~ 0.088), compact (P 01 -b ~ 2.45 h), highly inclined 
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(; ~ 87°), strongly relativistic, and displays eclipses of A and 
very large but systematic flux variations of B. Its proximity, 
relatively high measured flux density, and rapidly moving pul- 
sars (orbital velocities ^300 km s -1 ), make it an ideal target 
for Viss studies. In this paper we report measurements of the 
orbital modulation of V ISS from A at 820- and 1400-MHz us- 
ing data from the 100-m Green Bank Telescope (GBT). 

2. OBSERVATIONS AND DATA REDUCTION 

In 2003 December, our group was awarded 5x6hr Ex- 
ploratory Time tracks on 0737 with the GBT as part of 
the NRAO Rapid Science program. Three of the observa- 
tions discussed in this paper were mad e with the Berke ley- 
Caltech Pulsar Machine (BCPM; e.g. iKaspi et al1l2004 us- 
ing summed IFs and 72 /xs sampling; one at 1400 MHz us- 
ing 96 x 1 MHz channels and the other two at 820 MHz using 
96x0.5 MHz channels (Table [T). During one 820 MHz ob- 
servation, we also obtained data with the GBT Spectrometer 
SPIGOT Card, a new, flexible, and high performance pul- 
sar backend developed at Caltech and NRAO. It processes 
auto-correlations from the GBT Digital Spectrometer using 
two custom digital logic cards. These cards accumulate, sort, 
pack, and decimate the data before sending it to a PC for fur- 
ther packing and output into FITS files on a multi-TB RAID 
array. The SPIGOT currently handles one or two IFs with 
bandwidths of 50- or 800-MHz. We used a 50 MHz mode 
(#44) to process 1024 lags from summed IFs that were 8-bit 
sampled every 40.96 [is. In total, the SPIGOT observation 
generated 420 GB of data which we recorded onto removable 
hard disks. 

We folded these data modulo the predicted pul se period 
given the timing solution from iLvne et a"fl (120041) . These 
folds were used to determine a local timing solution for 
A and to create dynamic spectra from which to measure 
scintillation parameters. For the timing analysis, we mea- 
sured 300 times-of-arrival (TOAs) from each of the observa- 
tions (for individual integration times of ~ 1 min) by cross- 
correlating the folded profiles with a template profile from 
the SPIGOT data and then referencing the maximum of the 
cross-correlation with the time of the start of the observation 
as determined from the observatory clock. Our initial timing 
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FIG. 1.— Timing results from the 820 MHz GBT+SPIGOT observation 
listed in Table HI (Top) Residuals from a Keplerian orbit fit to the data. Sys- 
tematic deviations corresponding to the unmodeled Shapiro delay (shown as 
the line) are obvious in each orbit. (Bottom) Residuals from a fit using the 
"DD" timing model which includes the Shapiro delay for a 1 .25 Mg pulsar B 
at an orbital inclination i=88?4. The post-fit RMS residuals for the SPIGOT 
data alone are <7 fj,s. 

fits using TEMPO 8 showed approximately linear, frequency- 
dependent drifts in the TOAs of tens of fis during each of 
the 5-6 hr observa tions. Since A's p ulsations are significantly 
linearly polarized (Demorest et al. 2004), such apparent drifts 
may be caused by changes in the total-intensity pulse profile 
with time due to the rotation of the receiver feed with respect 
to the sky and gain differences in the two orthogonal IFs that 
were summed. In order to make high-precision measurements 
of local effects like the Shapiro delay, we fit for frequency 
jumps between the observations which effectively whitened 
the residuals to provide a post-fit RMS of ~ 1 1 [is. Using 
a mass for B of 1.25M© fLvne et a T| |i2004i ), we determine 
the Shapiro "shape" parameter to be sin; = 0.99962^ -{Jjj^f 
which corresponds to an inclination of 88.°4^j \. The errors 
for sin / and the other orbital parameters were determined us- 
ing a bootstrap analysis (e.g. lEfron"& Tibshirani 1986) which 
showe d that our results are consistent with those in lLvne et alJ 
(2004) at the 1-er level. Figure ^ shows residuals from the 
820MHz GBT+SPIGOT observation alone. 

Assuming a uniform Kolmogorov scattering medium, 
ICordes & Rickettl (119981) derived the velocity estimator 

Viss = 2.53 x 10 4 — ^kms" 1 , (1) 
v Atd 

where D is the distance to the source in kpc, Avd is the decor- 
relation bandwidth in MHz, Atd is the scintillation timescale 
in seconds, and v is the observing frequency in GHz. Note 
that unlike A?,/, which varies based on the relative velocities 
of the Earth, scintillation screen, and pulsar, Avd is a property 
of the screen itself and therefore does not vary significantly 
during an observation. The initial constant differs by factors 
of ~1— 3 depending on the structure function of the scatter- 
ing medium/screen as well as its location and extent. The key 
point is that Viss oc Atd' 1 , with an overall scaling that can be 
measured by fitting the known orbital velocity of the pulsar to 
the measured changes in Atd- 

In order to determine Atd and Avd for each of the obser- 
vations listed in Tabled we closely followed the process de- 
scribed in OBvS. Briefly, we created calibrated dynamic spec- 
tra (see fig. |3 by taking the on-pulse minus off-pulse flux 
during a duration 7j nt divided by the average level in each fre- 
quency channel during the full observation. We then com- 
puted the autocorrelation function (ACF) of non-overlapping 
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(and hence independent) blocks of the dynamic spectra of du- 
ration T^eas in order to measure Atd (defined as the l ie half- 
width of the central ACF peak in the time direction; ICordes! 
1986). We measured Avd (defined as the half-width at half- 
height of the central ACF peak in the frequency direction; 
lCordesl[l98 6) and its errors by fitting gaussians to the cen- 
tral ACF peaks from 15 sub-integrations and examining their 
statistics. We determined the relative errors for the Atd val- 
ues similarly, by examining the statistics for Atd fits made in 
6 and 3 frequency subbands for the SPIGOT and BCPM data 
respectively. Due to the small number of available subbands, 
we added additional errors of 5 and lOkm/s in quadrature to 
the measured errors for the SPIGOT and BCPM data respec- 
tively to guard against underestimates of Viss errors. 

We converted the measured Atd values into uncalibrated 
Viss measurements using eq. ^ ar, d calculated the true 
anomaly 6 at the center of eac h Viss measureme nt interval 
using the timing ephemeris from lLvne et alJ J2004I) . We per- 
formed a weighted least-squares fit of these data to eq. 7 in 
OBvS, Viss = n(v' 2 + v" 2 ) 1/2 , where v' and v" are the total ve- 
locities (i.e. orbital + systemic) parallel to and perpendicular 
to the orbital line of nodes in the plane of the sky respec- 
tively. In contrast to OBvS, we used the angle of periastron u> 
from the timing ephemeris and fit for four parameters; the sys- 
temic velocities parallel to (Vpi an e) and perpendicular to (Vp er p) 
the line of nodes, the inclination i, and the scaling parame- 
ter k. Fitting for k accounts for uncertainties in the nature of 
the scattering medium and the distance to the pulsar. Tabled 
lists the fitted parameters and their 95% statistical confidence 
limits as determined from projections of the x 2 space. The 
knowledge that we observe 0737 nearly edge-on (i ~ 90°) 
re moves the de generacy in possible inclination angles seen 
bv lLvnd (fT984) and OBvS and produces a single solution for 
each observation. The errors in the fitted parameters were cal- 
culated after scaling the Atd errors such that the reduced-^ 2 
of the fits equaled one. Combining the BCPM and SPIGOT 
results provides our best estimates of Vpi ane ~96.0±3.7kms _1 
and Vperp — 103.1±7.7kms _1 for a total systemic velocity in 
the plane of the sky of Viss ~140.9±6.2kms _1 . 

We note that Viss fits can result in i > 90°. In that case 
(as we have measured), the direction of V perp is either aligned 
with the orbital angular momentum vector if it is pointing 
slightly towards us, or anti-aligned if the angular momentum 
is pointing slightly away from us. Combining the timing- and 
scintillation-based measurements gives an inclination (in the 
more standard range 0°-90°) of i = 88.°7 ± 0.°9. 

2.1. Possible Biases to the Measured Viss 

While the "self-calibrating" nature of binary pulsar Viss 
measurements eliminates many of th e uncert ainties found in 
Viss studies of isolated pulsars (e.g. ICordesiri986l) . the mo- 
tion of the Earth with respect to the pulsar and the differential 
rotation of the Galaxy both contribute to the measured Viss- 

A flat, 220km s" 1 , Galactic rotation curve implies a ve- 
locity component in the plane of the sky of AV g j_ ~ 
14(D/0.6kpc)kms _1 due to differential Galactic rotation in 
the direction of 0737 Q = 245 .°2, b = -4.°5). However, as 
ICordes & Rickettl Jl998) point out, the ISM rotates along with 
the Galaxy, and we calculate that the effective transverse ve- 
locity is <lkms _1 for all reasonable distances to 0737. 

Earth's orbital motion projected onto the sky towards 0737 
contributes a doubly periodic velocity component AV©i_ to 
the measured Viss- For 0737, AV®j_ has minimum and max- 
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TABLE 1 

GBT Observations and V lS s Model Fits for J0737-3039A 





Start 


Freq. Res. 


Tabs 


Tim 


Tmeas 


Av d 


Viss 


yplane 


Vperp 


! 




X 2 /DOF 


(MHz) 


MJD 


(bands X MHz) 


(hr) 


(s) 


(s) 


(MHz) 


(km/s) 


(km/s) 


(km/s) 


(deg) 


K 


820 a 


52997.214 


1024x0.04883 


5.00 


10.0 


320 


0.096(4) 


139.0(4.5) 


95.1(3.8) 


101.4(8.0) 


91.0(1.1) 


0.631(11) 


0.64 


1400 b 


52984.195 


96x1.0 


6.12 


12.24 


612 


1.8(2) 


163(16) 


107(13) 


123(27) 


92.9(3.5) 


0.693(48) 


1.50 



NOTE. — " SPIGOT data. * BCPM data. All confidence intervals are 95% statistical. The number of degrees of freedom for the model fits are DOFg20 = 56-4 = 52 and DOF1400 = 
36 — 4 = 32. For these measurements, i can range from 0-180° (see i2*l . The errors on the fitted parameters were determined from projections of the x 2 space after adjusting the errors 
of the Viss measurements such that x 2 /DOF=l. 




FIG. 2. — Dynamic spectra (top) and Viss model fits (bottom) for the observations listed in TableQ The 820 MHz SPIGOT data is on the left and the 1400 MHz 
BCPM data is on the right. For the SPIGOT dynamic spectra, only the top 1/4 of the analyzed band is plotted. Values ±2-<r from the average level were set to 
black/white respectively in order to enhance the contrast. 



imum values of 23 and 30kms _I respectively. During these 
observations, AV®_l ~ 29kms _I . 

Since we do not know the orientation of the measured Viss 
on the plane of the sky, we cannot currently remove either of 
these biases from our data. However, measurement of an an- 
nual variation in Viss would allow us to remove AV^j_. Alter- 
natively, VLBA determination of the proper motion direction 
will allow us to subtract both biases directly. In addition, the 
combination of the VLBA angular proper motion [i with the 
transverse velocity V± ~ Viss will provide a unique geometric 
distance to the pulsar (D = V±/fi). 

3. DISCUSSION 

The relatively large Viss — 140kms _1 measured for 0737 
has important implications for the formation of the system. 
Here we assume that 0737 followed the standard formation 
scen ario for double neutron stars (DNSs) in the Galactic disk 
(e.g. |Bhattacharva & van den HeuveD ll991l) . In this picture, 
the first NS (A) forms in orbit about a massive (>8M Q ) stel- 
lar companion. The binary acquires some velocity due to 
impulsive mass loss in the supernova explosion and a pos- 
sible natal kick to the NS. At this stage, the system resem- 
bles a high-mass X-r ay binary (HMXB). Detailed studies 
(e.g.. lPfahl et al.l2002l) find expected HMXB systemic speeds 
<30kms _1 , because the NS kick is distributed over the large 



total mass of the binary. It is unlikely that the first supernova 
provided the high systemic velocity implied by the measured 
Viss- 

At the end of the HMXB phase of DNS formation, the 
stellar companion evolves, fills its Roche lobe, and trans- 
fers matter to the NS. Because of the extreme mass ra- 
tio of the binary, the transfer is dynamically unstable, and 
leads to a common-envelope phase. If coalescence is 
avoided, the NS emerges in a tight orbit with the hydrogen- 
exhausted core of its companion. The stellar core contin- 
ues its nuclear evolution, and may stably transfer matter to 
the NS via Roche-lobe overflow. This is the favored mech- 
anism for spinning up and reducing the d ipole magnetic 
field strength of A (Dewi & van den Heuvel 2004). Subse- 
quently, the stellar core explodes and produces the second 
NS (B). In 0737, the pre-supe rnova mass of the stellar core 
was likely M r ~ 2.0-3.5 M r, (iDewi & van den HeuvellkOO^ 
IWillems & Kalo gera 2004|). The combination of supernova 
mass loss and a kick to B could easily have given 0737 a 
speed of >140kms _1 . We now use the interestingly large 
component (Vperp ~103kms _I ) of the systemic velocity per- 
pendicular to the orbital plane to constrain the magnitude of 
the kick imparted to B. 

The masses of A and B are M A ~ 1.34M and M B ~ 
1.25 Mq respectively, giving a total binary mass of M t = 
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M A + M B ~ 2.59M Q l lLvne et al.H2004 . We assume that the 
pre-supernova orbit was circular, and denote by R, V, and 
H = R x V the position, velocity, and orbital angular mo- 
mentum of the pre-collapse progenitor of B relative to A. If 
B receives an impulsive natal kick velocity of V*, the post- 
supernova relative orbital velocity and angular momentum are 
given, respectively, by V = V + V* and H' = H+R x V*. The 
velocity acquired by the binary center of mass after the super- 
nova is 

M A (M C -M B ) M B 

Vcm = -TT 7iF J\ V+ TT V *- <2> 
M, (M C +M A ) M, 

If Vk = 0, then Vcm is in the orbital plane. Therefore, the 
large observed perpendicular component of Viss, Vperp, ma Y 
require B to have been born with a significant kick velocity. 
If we assume that Vcm, and thus the observed V P erp, is entirely 
due to the second supernova explosion, it is straightfor ward to 
show that Vp erp is related to V* by (see also Wex et al. 2000) 



H' V 



CM : 



: H'V 



'H-V* 1 + 



Ma 



(3) 



For specified masses and orbital parameters before and after 
the supernova, Vk is minimized when V* is parallel to H (i.e., 
perpendicular to the pre-supernova orbital plane). The mini- 
mum kick speed is then 



V k>B 



■Vr 



perp 



~H 



i + 



Ma 

M r 



Note that 
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M, 



Mr+Mi 



-d-e' 2 ) 



1/2 



(4) 



(5) 



and <?' are the immediate post-supernova semi- 
major axis and eccentricity, respectively, which have since 
evolved to their observe d values under the action of gravi- 
tation al wave emission. iDewi & van den He uvel (2004|) and 
IWillems & Kalogera! (12004 each find that e' < 0.14. so that 
a' I a cannot differ substantially from unity. For our purposes, 
it is sufficient to assume that (a'/a)(l-e' 2 ) ~ 1. Finally, we 
obtain 



^ I Mr 



1/2 



1 + - 



Ma 

Mr 



1/2 



(6) 



For the range of pre-collapse core masses given above, 



and V m 



dOOkms , we find that Vk, 



100-150kms" 



where the upper limit corresponds to M c = 2M Q . If the 
pre-supernova core mass takes its smallest possible value of 
M c = M B , then Vk^ \ « — 21 3 km s" 1 . 

Piran & Shaviv (2004) point out that at a distance of 
0.6 kpc, 0737 is only ~50pc from the Galactic plane. They 
suggest that in order to find 0737 so close to the plane, Vcm 
must be very small (<15kms _1 or <150kms _1 at 68% or95% 
confidence respectively), and that the supernova mass loss and 



kick to B must likewise have been very small. The measure- 
ment of a large Viss implies that (1) either the velocity direc- 
tion is largely in the Galactic plane or we happen to see it as it 
passes through the plane, and (2) moderate pre-collapse core 
masses and NS kicks are permitted if not required. 

VLBA observations will soon determine the mag- 
nitude and direction of 0737's proper motion [i ~ 
0.05"yr I (Viss/140kms- I )(D/0.6kpc)- 1 . For long-term tim- 
ing observations, the large velocity in the plane of the sky 
V i ^ Viss will cause an apparent acceleration dShklovskiil 
1970) that will bias the measured spin- and orbital-period 
derivatives. For 0737, P/P = V|/(Dc) = V±n/c ~ 3.5 x 
10~ 18 s" 1 , which amounts to ^5% of A's measured spindown 
and ~2.5% of the predicted orbital decay rate due to gravita- 
tional wave emission. Future bias-corrected measurements of 
Viss should determine V± to a few percent and VLBA observa- 
tions will likely measure [i to even better precision. Therefore, 
the "Shklovskii effect" should be measurable to <5% preci- 
sion and consequently, its error will contaminate the measure- 
ment of the orbital period-derivative P 01 -b due to the emission 
of gravitational wave emission by <0.2%. 

4. CONCLUSIONS 

We have measured the systemic velocity and inclina- 
tion of 0737 using the "self-calibrating" method of binary 
scintillation velocity measurements. The inferred high ve- 
locity Viss— 140kms _I strongly suggests that a substantial 
(>100kms _1 ) kick was imparted to pulsar B at its birth. The 
large Viss will impact long-term timing of the system and will 
allow a precise geometric distance measurement when com- 
bined with a VLBA proper motion. 

Future scintillation observations using the high frequency 
resolution and wide bandwidth of the GBT+SPIGOT will 
allow substantially improved measurements of Viss at 1400 
and even 2200 MHz where we have already detected the or- 
bital modulation of A in BCPM data. These measurements 
should also allow the removal of the contaminating effects of 
the Earth's motion and the differential Galactic rotation. Fi- 
nally, we have already detected scintillation from B during the 
bright portions of its orbit, but future observations may allow 
"snapshot" calibrations of Viss based on measurements of Af</ 
from A and B at an instant in time given our knowledge of the 
relative orbital velocities of the pulsars. Such measurements 
would demand far less telescope time. 
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